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We construct a quantum algorithm that performs function-dependent phase transform and re- 
quires no initialization of an ancillary register. The algorithm recovers the initial state of an ancillary 
register regardless of whether its state is pure or mixed. Thus we can use any qubits as an ancillary 
register even though they are entangled with others and are occupied by other computational pro- 
cess. We also show that our algorithm is optimal in the sense of the number of function evaluations. 

PACS numbers: 03.67.Lx, 03.65.Bz, 89.70+c 

Quantum computation is based on three quantum phenomena: superposition of states, quantum interference, 
and quantum entanglement. These effects enable exponential speedups in the solutions of certain problems and 
allow one to transgress some boundaries of classical computational complexity theory (l - 1 1 1 . Most known quantum 
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algorithms (l|||,^,^,^T[^2) rely on conditional phase transform the realization of which is accomplished by the quantum 
Fourier transform or the Walsh-Hadamard operator together with the unitary operator evaluating a given function 
or quantum oracle. In general, conditional phase transform can be described by the operation ^2 x —q ol x \x) \— > 
12x=o ex P[2' n "if( x )/M}a x \x) for a function / : 1^ — ► Zjf, which we call / -dependent phase transform. The resulting 
interference pattern facilitates determining global property of the underlying function. 

• We need a quantum circuit evaluating a function to perform function-dependent phase transform and unitary 
| evolution of quantum computational process requires an ancillary register from which we have to extract the desired 

• relative phases conditioned on the given function. All previous quantum algorithms resort to initialization of the 
ancillary register before the computation. We may ask a question: Is it possible to perform function-dependent 
phase transform without initializing and deforming the state of the ancillary register? If it were possible, energy 
dissipation caused by initialization process would be avoidable and any register that could contain useful information 

■ to be preserved could temporarily be used as an ancillary register. In this work we construct a quantum algorithm 
that implements function-dependent phase transform without initializing an ancillary register. Furthermore, the 
application of the constructed algorithm retrieves the initial state of the ancillary register. Thus the ancillary register 

O || can consist of any qubits collected from any other registers even though they are being used in other computation 

■ which can proceed after carrying out their auxiliary duty in function-dependent phase transform. We show that to 
realize function-dependent phase transform at least two operations dependent on the given function are necessary. 



Thus the presented algorithm is optimal in the sense that it involves only two function evaluations. Of course, if any 
kind of initialization is involved, one function evaluation is sufficient. 

The /-dependent phase transform TZkj ■ \%) l— * w m \ x ) P^ a y s an important role in quantum algorithms where 
. , u>m — exp(27ri/M) is a primitive M-th root of unity and k € Zm may be chosen appropriately depending on the 
given problems. For simplicity, we assume that N and M are powers of 2, that is, N — 2" and M = 2 m for some 
nonnegative integers n and m. In order for the information on the given function to be encoded in the phases it is 
necessary to evaluate the given function on quantum computer. On quantum computer the evaluation of a function 
is performed by the unitary operation Uf : \x) ® \y) i— > \x) ® \y + f{x)) for x £ 1>n and y € Zm- The first n-qubit 
register we call the control register contains the states we wish to interfere. The second m-qubit register called the 
function or ancillary register is used to draw relative phase changes in the first register. The superposition principle 
of quantum mechanics allows us to prepare the computer in a coherent superposition of input states and to compute 
exponentially many values of / in superposition with a single application of Uf. This phenomenon is the basis for 
quantum parallelism which leads to a completely new model of computation. In view of the second register the 
function evaluation adopts a translation operator T z : \y) t— > \y + z) where z is dependent on the state of the first 
register. The operator Uf can be described in terms of the translation operator on the second register; 

u f ■ ax v 1^ ® ly) ^ H a *v \ x ) ® T /(z) \v) ■ W 

x.y x,y 

We note that to implement IZkj at least two registers are necessary due to Uf. By adding an ancillary register the 
effect of the li-kj on the control register can be viewed as phase changes in the ancillary register dependent on the 
states of the control register. To be more specific, we define a unitary operator J^,z '■ \y) l— * w m \y) f° r y € I^m- Then 
T^k f ® I can explicitly be written by 
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TZ k j ®1 a xy \x) (g> \y) h-> ^ a XJ/ \x) ® J k j( x ) \y) ■ (2) 

x.y x,y 

We note that J k:Z has one eigenvalue wff and the corresponding eigenspace is the whole Hilbert space. Due to the 
expressions ([!]) and (|^) we can concentrate on the operations of the ancillary register. 

Especially when / is the identity map T, TZk,i maps \y) to \y), in which the phase-encoded information depends 
on its state, and can be obtained by T^TuT where T is the quantum Fourier transform. A quantum algorithm 
to implement Jk,z can be realized using 7?.fc,r and 7^. We prepare an arbitrary m-qubit register whose state is 
I*) = Tl,y=v a v \v) an d proceed the following algorithm: (i) Apply T z . (ii) Apply TZ k j. (iii) Apply Tj = T_ z . (iv) 
Apply 1Z\. j = T^T k T . Then the state evolves as follows: 
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y=o 

«>% m ■ (3) 

Therefore we get 1Z\ T T^lZ k .zT z Y$>) = uj^f for an arbitrary namely, we have 

j h , z = J^t = nl x T^n k .iT z . (4) 

Theorem 1. There exists a quantum algorithm to implement J k ^ z using two T± z . 

The algorithm for J k , z is not unique. All cyclic permutations of the steps are identical. For example, we can start 
at Step (ii), perform successive steps, and end at Step (i). In fact, if we use the notation [A, B] = ABA~ 1 B~ 1 then 
by Eq. (^) we have 

Jk,z = [nl I ^] = [r z ,nl x ] 

= [Kk,i,T z } = [T2,n ktI } (5) 

with its inverse J-k, z = [TZk,i,T z ] = [%,T^k,i] = \R-\,i,T z ] = \T),v} kT \. Furthermore, noting that S k ,i = FT~ k F 
maps \y) to \ —y) one can easily check that 

Jk lZ = Sk,xT z Sk,iT z , (6) 



hence the overall unitary operation performed by the algorithm can be written by V2^± z Vi = for some unitary 
operators Vi and V2. Since ViV^T^z — ^^1, it is enough to consider a unitary operator V such that V7± 



which also offers another implementation. We remark that Sj, x = S k ,x- Therefore there are many methods to 
implement J k . z . However, the number of T± z in each implementation is always equal to two and cannot be reduced. 
Let us suppose that there exists a quantum algorithm implementing J k ^ z where the only way to implement the 

given information on z is through T± z . Then the dependence of z requires at least one T± z at a certain step and 

kz e - 

Af J 

Since V = <^mT± z -, V depends on z. Thus in some another step of the algorithm we have to use information on z 
once more and so the overall procedure includes at least two translations by ±z. This observation will later be used 
in showing that our algorithm for 7t k j 1S optimal. 

We now turn to the /-dependent phase transform TZ k j- Wc let |$) = J2x=o ax \ x ) 

and |*) = E;=o Pv\v) be the 

respective states of the control and the ancillary registers. It is noted that no initialization is involved during the 
preparation of the registers. By inspecting Eqs. (|l|) and (Q) the algorithm (||) for J ktZ leads to an algorithm for TZ k j: 
(i) Apply Uf. (ii) Apply I ® H kt x. (iii) Apply Ut = U-f. (iv) Apply X ® TZ\ T =1 ® 1Z- ky x. This procedure makes 
the state of the registers evolve as follows: 
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JV-1 M-l 

|$}®|*) X) E a -A/I*>® ly + /(«)) 



a;=0 j/=0 
JV-1 M-l 

£ ^ a:c /3rf +/(x)) |x)®|y + /(x)) 

a;=0 y=0 
, t N-l M-l 



x=0 y=0 

kf(x) 

x=0 



Now we discard the ancillary register. Then we obtain the /-dependent phase transform TZkj ■ \x) i— > \ x )- 

Theorem 2. There exists a quantum algorithm that implements function-dependent phase transform using two 
evaluations of a given function such that the ancillary register preserves its initial state. 

Since the /-dependent phase transform TZkj can be written in terms of Jk,z as in Eq. (§) and by Eq. (§) there are 
many methods to realize SJk,z-> we can conclude that the algorithm for TZkj is also not unique. 

In the procedure (^) we have assumed that an ancillary register is in a pure state. However, this is not an essential 
requirement. In fact, any mixed state is allowed. To be more precise, let A be a quantum system to be used as 
an ancillary register and its state be described by the density operator p A . Then there exists a reference system 
R such that the compound system AR is in pure entangled state |\I>- A - R ) that gives rise to the given reduced state 
p A = Tir(p ar ) where p AR — \ty AR } (\ft AR \ is called purification of p A . Using the Schmidt decomposition we can 
rewrite | xf" 4 ^^) as X^lo* a v \v A ) ® l^^) - ^ e n °t e that the states |\&*)'s may not form the standard basis for the 
subsystem R but just an orthonormal basis while the states |y A )'s form the standard basis for the subsystem A. Now 
applying the above algorithm to |$) (g> |<I'" 4fi: ) one can see that the final state becomes (TZkj |$)) <8> V? ). Thus our 
algorithm works whether the state of the ancillary register is pure or mixed. This implies that we can compose an 
ancillary register of any m qubits which are collected out of any other registers even though they are still being used 
in other computational process and are possibly entangled with other qubits. The presented algorithm ((7J) recovers 
the initial state of the joint system AR after extracting the desired relative phase changes. Thus the qubits in the 
temporarily composed register can be restored to their positions to continue the suspended computation. 

Our algorithm requires two evaluations of /, i.e., Uf and U\ [or two Uf when Eq. (^) is applied]. This is because 
we employ no initialization. We know that Uf causes translations in the ancillary register by Eq. (|l|) and that any 
quantum algorithm for J7fe i2 adopts at least two translations. Thus we have the following theorem. 

Theorem 3. Any quantum algorithm that implements function-dependent phase transform without initialization 
requires at least two evaluations of a given function. 

On the other hand, if the ancillary register is initializable only one evaluation of / is sufficient. Indeed it is clear 
that J-\—k) is an eigenvector of T z with the corresponding eigenvalue w^f . If we let |\&) = TT-k |0), then Uf maps 
I a;) ® I*) to w'f 1 ' |x> ® I*). The special case for k = 1 was studied in Jl|@. 

Initialization in general sense is a process to transform the state of a quantum system to a definite pure state, 
which can later be rotated to |0) as usual or any other desired state by frame change. When we are to initialize the 
subsystem A we cannot avoid corrupting the correlation between the subsystems A and R, which can be measured 
by the quantum mutual entropy S(A : R). If the subsystem A is entangled with the reference system R, that is, 
S(A : R) ^ 0, then even when p A is known we cannot initialize the subsystem A by local unitary operations on A. We 
note that the quantum mutual entropy is invariant under local unitary operations of product form for each subsystem. 
If the bipartite systems A and R are separable, that is, S(A : R) = 0, then a certain frame change on the subsystem 
A effects on initialization of the ancillary register without knowing the total state p AR . In this sense we say that the 
subsystem A is nondestructively initializable when p A is pure and known and initialization by a local frame change on 
the subsystem A is called nondestructive initialization. If nondestructive initialization is adopted then the ancillary 
register regains its early state and one evaluation of a function is sufficient for function-dependent phase transform. 

Let us consider a more general function / : Z^v — > [0, 1) Cl. We define m-bit approximation / : Zjv — > %m of / 
by f( x ) — y^i—i a-i2' m ~ l G for x £ Zjv where (0.aia2 . . . a m )2 = X)™ 1 a »2 _I is an m-bit binary expansion of f(x) 
for a, G Z2. Then the approximate f -dependent phase transform lZ k ? approximates the operation \x) 1— > e 2 ™ 1 ?^ 



||l3| . [l4| . This approximate /-dependent phase transform can be applied to the conditional 7-phase transform and the 
/5-phase diffusion transform in [fL2|. 
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All known quantum algorithms are based on the effect of function-dependent phase transform and the quantum 
Fourier transform (or the Walsh-Hadamard operator) . The quantum Fourier transform enables one to find the period 
of a function in polynomial time and plays an essential role in Shor's quantum polynomial-time algorithms |(| for the 
integer factoring and the discrete logarithm problems which are known to be intractable on classical computer. It 
also enables us to construct function-dependent phase transform without initializing an ancillary register, which can 
immediately be applied to most known quantum algorithms. 

Deutsch and Jozsa |l],|2| presented a simple promise problem to determine whether a Boolean function / : — > Z2 
is either constant or balanced and showed that it can be solved efficiently without error on quantum computer while 
it requires exhaustive search to solve deterministically without error in a classical setting. The key of their algorithm 
is the 7r-rotation of phases controlled by the query result of quantum oracle. In this problem M = 2 and k = 1. 
Then T becomes the Walsh-Hadamard operator ^ = 775 ( 1 _i) an d T\ = T-i becomes the Pauli spin operator 

o~ x = ^ 1 0) wn ich represents a bit-flip. The operator TZi,x is a phase-flip operator a z = -1)' Thus the overall 
scheme for 1Zij is (J <E> cr z )Uj (I <£> cr z )Uf |lS{ ]. Grover constructed a quantum algorithm that can find a particular 
item in expected time 0{\fN) when an unstructured list of N items are given His algorithm relies on the 

conditional phase transform Sj : \x) i— > (—1)^^' \x) and the diffusion transform D — WSqW^ where / is the Boolean 
function computed by an oracle, So = Sf , and fo(x) — 6q x - In this case Sf — IZij with M — 2. Brassard and 
H0yer |[l| combined Simon's ZQP algorithm || and Grover's quantum search algorithm and showed that Simon's 
problem can be solved on a quantum computer in worst-case polynomial time and thus is in QP class. Thus their QP 
algorithm mainly depends on the conditional phase transform. Chi and Kim |i~2[ generalized Grover's algorithm and 
showed that a quantum computer can search a database by a single query when the number of solutions is equal to or 
more than a quarter. Their algorithm makes use of the conditional 7-phase transform Sf n 

; \ x ) I — > e <7/(*) \ x ) 

and the 

/3-phase diffusion transform Dp — WjSj^Wj where Wi is any unitary transformation satisfying Wi \l) — J2x=o \ x ) 
and Si t p = Sf lt p with fi(x) = §i x - In this case we can use approximate function-dependent phase transform. 

All conditional phase transforms fall into the category of function-dependent phase transform. Therefore our 
algorithm for function-dependent phase transform is directly applicable to most quantum algorithms. We note that 
for general positive integers N and M the approximate Fourier transform in [jl6f can be used in our algorithm. 

In summary, we generalized conditional phase transform to function-dependent phase transform and presented a 
quantum algorithm that performs function-dependent phase transform and does not require any kind of initialization 
of an ancillary register. Our algorithm recovers the initial state of the ancillary register. Thus we can compose 
an ancillary register of any qubits regardless of whether they are entangled with others or being used in another 
computational process. Our algorithm employs two evaluations of a given function and is optimal in that any 
quantum algorithm that implements function-dependent phase transform without initialization requires at least two 
evaluations of a given function. 
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